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ABSTRACT 

We have discovered that the spectrum of the well-known dwarf nova EM Cyg 
is contaminated by light from a K2-5V star (in addition to the K-type mass donor 
star). The K2-5V star contributes approximately 16 per cent of the light from the 
system and if not taken into account has a considerable effect upon radial velocity 
measurements of the mass donor star. We obtain a new radial velocity amplitude 
for the mass donor star of K2 — 202 ± 3kms _1 , which compares with the value of 
K2 — 135 ± 3kms _1 obtained in Stover, Robinson & Nather's classic 1981 study 
of EM Cyg. The revised value of the amplitude combined with a measurement of 
rotational broadening of the mass donor v smi = 140 ± Gkms" 1 , leads to a new 
mass ratio of q = M2/M1 — 0.88 ± 0.05. This solves a long standing problem with 
EM Cyg because Stover et al.'s measurements indicated a mass ratio q > 1, a value 
which should have led to dynamically unstable mass transfer for the secondary mass 
deduced by Stover et al. The revised value of the mass ratio combined with the orbital 
inclination i = 67±2° leads to masses of 0.99±0.12M Q and 1.12±O.O8M for the mass 
donor and white dwarf respectively. The mass donor is evolved, since it has a later 
spectral type (K3) than its mass would imply. 

We discuss whether the K star could be physically associated with EM Cyg or 
not, and present the results of the spectroscopic study. 

Key words: binaries: spectroscopic - novae, cataclysmic variables - accretion, ac- 
cretion discs - techniques: radial velocities 



1 INTRODUCTION 

EM Cygni is one of the most familiar examples of the 
sub-group of cataclysmic binary stars called dwarf novae 
(DN). They consist of a white dwarf star accreting material 
through the inner Lagrangian point from a late spectral- 
type mass donor. The material spirals down towards the 
white dwarf and forms an accretion disc around it. EM Cyg 
belongs to the sub-class of dwarf novae named after the pro- 
totype Z Cam, in which the peculiar characteristics of the 
light curve were first observed. They experience - in addi- 
tion to the normal outburst/quiescent states of DN - peri- 
ods of time when they appear to be 'stuck' in a semi-high 
state, making them appear visually approximately 0.7 mag 
brighter than at minimum. 

Additionally, EM Cyg is one of only a handful of cata- 
clysmic variable (CV) stars which show an eclipse of the disc 
every orbit and spectral lines from both stellar components 



of the binary. This makes EM Cyg special amongst CVs, as 
it means that accurate binary parameters can be deduced 
from spectroscopic and photometric observations. In a clas- 
sic study Stover et al. (1981) measured its absorption line 
radial velocity amplitude together with that from the emis- 
sion lines and obtained a mass ratio q = M2/M1 — 1.26; 
this value confirmed an earlier study of Robinson (1974), 
who calculated that q — 1.29. A mass ratio greater than 
one is remarkable because if mass transfer is from the more 
massive star to its companion, the orbital separation must 
shrink in order to conserve angular momentum. Unless the 
donor star shrinks faster than the Roche lobe, mass transfer 
is unstable. The instability occurs on a dynamical time scale 
if the adiabatic response of the star fails to keep it within 
its Roche lobe. At higher donor masses (M2 > O.8M0) how- 
ever, the star shrinks drastically when it loses material at 
a high rate and so the mass transfer becomes stable on a 
dynamical time scale but it is still unstable thermally (e.g. 
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de Kool, 1992, Politano, 1996). The mass ratio (q=1.26) and 
the donor mass (Ma = 0.76M Q ) of Stover et al. (1981) lead 
to dynamical instability according to the criteria given by 
Politano, 1996: EM Cyg should not exist. 

Despite this unique status, there have been no more re- 
cent studies of EM Cyg. In part, this is probably because 
measurements of the radial velocities of the white dwarfs 
in cataclysmic variables are normally based upon the emis- 
sion lines which usually come from the accretion disc which 
surrounds the white dwarf, and these are well known to be 
unreliable (Stover et al., 1981). It would perhaps not be sur- 
prising if this problem affected Stover et al.'s (1981) study 
of EM Cyg. However, EM Cyg is one of the longer period 
dwarf novae (P = 6.98 h), for which the distortions of the 
emission lines are less significant because of the relatively 
high amplitude of the white dwarf, and in Stover et al.'s 
study, the emission and absorption line amplitudes are only 
20° from being in anti-phase, by no means a large distortion 
for these stars. 

In this paper we present new spectra of EM Cyg which 
show clear evidence for contamination from a third star, 
which happens to have a very similar spectral type to the 
mass donor star. Accounting for the presence of this star in- 
creases the radial velocity semi-amplitude of the mass-donor 
star considerably, and causes the mass ratio (q = M2/M1) 
to become less than one, at the same time leading to a larger 
donor mass. Thus the mass transfer in EM Cyg is therefore 
stable as expected from theory. 



2 OBSERVATIONS 

On 22 June 1997 we took 102 spectra with the Intermediate 
Dispersion Spectrograph (IDS) on the 2.5 m Isaac Newton 
Telescope (INT) on the island of La Palma in the Canary 
Islands. Our dispersion was 0.4 A per pixel and the spectra 
covered 6230 to 6650 A. The resulting spectra have a reso- 
lution of 0.8 A (FWHM). Conditions were clear throughout 
the night and the seeing was approximately 1 arcsec. 

Each target exposure was of 200 s duration. The 0.8 
arc sec wide spectrograph slit was oriented at a position 
angle of 255° to capture the spectrum of a star 22 arcsec 
from EM Cyg in order to calibrate slit losses. Arc calibra- 
tion spectra were taken every 35 minutes or so. The arc 
spectra were fitted with fourth-order polynomials. The rms 
scatter of these fits was of the order 0.003A. Arc spectra 
were linearly interpolated in time to provide the wavelength 
calibration for EM Cyg. The CV appeared to be approach- 
ing its low state after a standstill when the spectra were 
taken (Mattei, 1999). Fig. [j] shows the light curve behaviour 
of EM Cyg during the period April 1997 to February 1998. 

The CCD frames were bias-subtracted and then cor- 
rected for any pixel-to-pixel sensitivity variations using ex- 
posures of a tungsten lamp. Variations in the illumination 
of the slit were corrected using observations of the twilight 
sky. We extracted the spectra with optimised weights (Home 
1986). Spectra of radial- velocity standard K and M dwarf 
stars were also obtained. In all, 25 standard stars were ob- 
served, selected from the lists of Marcy et al. (1987), Beavers 
and Eitter (1986) and Duquennoy, Mayor and Halbwachs 
(1991). 
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Figure 1. The light curve of EM Cyg during 1997. The arrow 
marks the date of the observations. 
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Figure 2. The mean spectrum of EM Cyg observed on 22 Jun 
1997. 



3 RESULTS 

The mean spectrum of EM Cyg is presented in Fig. |^ (lower 
spectrum). The upper spectrum in Fig. ^ shows the aver- 
age of the 102 spectra after subtraction of the contribution 
to the absorption lines from the co nta minating K-type star 
(hereafter 'third star'), see section 3.1, and after each spec- 
trum has been shifted into the rest frame of the mass donor 
in EM Cyg in order to remove the effect of orbital smear- 
ing. Absorption features characteristic of a K-type star are 
visible in the lower spectrum. They are sharp, which is un- 
expected because the lines from the mass donor should be 
broadened by rotation, and by orbital motion. These lines 
are from the third star. It can be seen from this plot that 
the absorption lines in the upper spectrum (which represent 
those from the mass donor of EM Cyg) are indeed much 
wider than those in the lower spectrum, which is as expected 
due to the broadening effects previously noted. 

This is confirmed from the trailed spectrum presented 
in Fig. |. 

The procedure followed to get from the lower spectrum 
to the upper spectrum in Fig. H was fairly complex, and so 
we start with a summary of the main points. 



(i) The continua of the individual EM Cyg spectra were 
fitted and then subtracted. The resulting spectra were then 
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Figure 3. The trailed spectrum of EM Cyg scaled to show the absorption features. The contaminating K star is visible in narrow lines 
that run straight up the plot. The donor star is scarcely visible, but sinusoidal behaviour can be spotted in the feature at 6495 A. Ha is 
not included on this plot (see Fig. fill). 



re-binned onto a uniform velocity scale. The radial-velocity 
standard stars also had their continua subtracted, and were 
then binned onto the same velocity scale. 

(ii) The average EM Cyg spectrum was cross-correlated 
(using the method of Tonry & Davis, 1979) with each dif- 
ferent radial-velocity standard star. This was done in order 
to obtain a list of the relative shifts between the two sets 
of absorption lines. These shifts were then applied to the 
standards, so that the absorption lines were congruous with 
those of the EM Cyg spectrum. 

(iii) The standard star spectra were then scaled and sub- 
tracted from the EM Cyg spectra in order to remove the 
absorption lines of the third star. We refer to this process 
as "optimal subtraction". 

(iv) The spectral type of the standard star which gave the 
lowest x 2 on subtraction was adopted as the spectral type 
of the third star (e.g. see Fig. |^) . 

(v) The individual, contaminant-free EM Cyg spectra 
were then cross-correlated again with the radial-velocity 
standards. This time the standards were artificially broad- 
ened t o m atch the rapid rotation of the mass donor (see 
section 3.2). The resulting radial velocities were then fitted 
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with circular orbits in order to determine the radial velocity 
semi- amplitude of the mass donor, and the systemic velocity 
of the EM Cyg. 

(vi) The fit which minimised \ 2 was then shifted out of 
the individual EM Cyg spectra in order to correct for the 
orbital motion of the mass donor. The resulting spectra were 
averaged to obtain a spectrum of EM Cyg in the rest frame 
of the mass donor (Fig. ^, upper spectrum) . 

(vii) In order to derive the best-fit spectral type, and 



match the v sini of the mass donor steps (i) to (iv) were 
repeated including the step of artificially broadening the 
radial- velocity standards (as in step (v)). See section 
for details. 



3.1 Analysis 

The presence of the donor star is hard to see in the raw data 
(Fig. ^) , but becomes clear after cross-correlation with a K- 
star radial- velocity standard (Fig. |^). The cross-correlation 
technique used is similar to that used for SS Cygni (Stover 
et al., 1980). The spectra had their continua fitted and then 
subtracted and were then transformed onto a logarithmic 
scale so that the cross-correlation technique could be used 
to measure the velocity shifts due to Doppler effects. In order 
to effectively remove the contribution to the absorption lines 
from the third star, the relative shifts between the velocity- 
standard stars and the average EM Cyg spectrum were de- 
termined, and then applied to the velocity standards. Then 
the standards were subtracted from the EM Cyg spectra to 
find the optimum frac tion at which the standard should be 
scaled; see section 3.3 for a more detailed description of the 



procedures followed. This routine also provides a formal er- 
ror on the fractional contribution from the third star, and 
the standard star which gives the minimum \ 2 indicates its 
optimum spectral type. The standards were then multiplied 
by this constant, and subtracted from the EM Cyg spectra. 
The optimum fraction obtained is 0.160±0.002, with a K3V 
standard star. 

The spectral type of the third star could be initially de- 
duced from the appearance of the absorption line spectrum 



© 1999 RAS, MNRAS 000, 000-000 



4 R.C. North, T.R. Marsh, C.K.J. Moran, U. Kolb, R.C. Smith, R. Stehle 




-500 

Velocity (kms -1 ) 



500-500 

Velocity (kms -1 ) 



500 



Figure 4. Trailed spectra of the cross-correlations with a K-type standard a) before subtraction of a K3V standard, left, and b) after 
the subtraction, right. As can be seen from a), the contaminating spectrum produces another set of correlation peaks, visible here as 
a straight line going up the trailed spectrum centred at -20 kms -1 . The correlation peaks measuring the radial velocities of the mass 
donor around the orbit are also visible in both plots. No artificial broadening has been applied to the standard in plot a), but in plot b) 
the standard has been artificially broadened to ■usini = 140 kms -1 



between 6350A and 6540A (see Fig. ||). The criteria devel- 
oped by Casares & Charles (1993) to determine the spectral 
type of the mass donor in V404 Cyg at wavelengths around 
Ha were also used to classify the contaminant. In particu- 
lar, the lines at CalAA6439.1, and FelAA6400.0 + 6400.3, 
and the blends at Cal AA6449.8 + 6450.2 and Ca I A6462.6+ 
FelA6462.7, whose relative depths are useful indicators of 
spectral type (Casares & Charles, 1993) were used. How- 
ever, as can be seen from Fig. |B| the relative depth of the 
blends at wavelengths 6400 A and 6440 A are also sensi- 
tive to variations in T e g. The blend at 6400 A dominates 
at spectral types around G8V but then diminishes relative 
to the blend at 6440 A becoming equal at around K3, un- 
til it is hardly apparent at spectral types around M4V. In 
EM Cygni these blends appear to have equal strengths rel- 
ative to each other, indicating a spectral type around K3V. 
The optimal subtraction process described earlier in this sec- 
tion was carried out for each individual EM Cyg spectrum, 
and a different constant computed for each. Then the stan- 
dards were scaled using these constants, and subtracted from 
the EM Cyg spectra. The spectral type of the third star was 
best fitted by standards between types K2V and K5V. The 
standard which minimised the value of \ 2 had a spectral 
type of K3V. 

The third star subtraction process was then carried out 
'in reverse' to check for possible systematic errors caused 
by this spectral typing procedure. This means that the ab- 
sorption lines from the mass donor in EM Cyg (which are 
broadened by rotational and orbital effects) were optimally 
subt ract ed first (using the same general method as in sec- 
and then the optimum spectral type of the third 



tion 



3.2) 



star was determined. In this case the optimum spectral type 
for the third star came out slightly earlier at K2V, indicat- 
ing that we cannot constrain the spectral type any more 
precisely than between K2V and K5V with this method. 



Finally the sequence of cross-correlation, shifting, and 
optimal subtraction was repeated using the best-fit spectral- 
type standard to the contaminating absorption lines (K3 V). 
The resulting spectra were then used (see section 3.2) to 
find out the spectral type of the mass donor in EM Cyg. The 
relative shift between the standard star absorption lines and 
those of EM Cyg's average spectrum gave a measure of the 
systemic velocity of the third star. This turned out to be 
7 = —20 ± 3 kms -1 . This is similar to the systemic velocity 
derived for the mass donor in EM Cyg by Stover et al. (1981) 
from the absorption lines of the spectrum in the H/3 region 
(— 23±6 km s -1 ), and is also similar to the value determined 
in section 3.2 (—25 ± 2kms -1 ). The agreement in systemic 



velocities of the mass donor and the third star is a hint that 



they are physically associated (see section 4.2 ). 

The next section details the method used to obtain the 
radial velocity semi-amplitude of the mass donor (which in- 
corporates the cross-correlation technique followed by Stover 
et al. (1980)). 



3.2 Radial velocity of the Donor star 

Once the individual spectra of EM Cyg had been corrected 
for the presence of the third star, the radial velocity semi- 
amplitude of the actual mass donor was calculated. The 
corrected EM Cyg spectra were cross-correlated (using the 
method of Tonry & Davis, 1979) with the radial-velocity 
standards each of which had been artificially broadened. The 
resulting radial velocity curves were then fitted with circular 
orbit fits of the form 



V{t) = 7 + Ksva 



Then each EM Cyg spectrum was shifted to correct for the 
orbital motion of the donor star, and the results averaged. 
This gave the average spectrum of EM Cyg in the rest frame 
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Figure 5. The average spectrum of EM Cyg scaled to show the 
absorption lines in the wavelength region of 6350 - 6540A , placed 
in an apparently suitable position for the third star between spec- 
tral types K3V and K5V. The broadened absorption lines from 
the mass donor underlie those from the third star. The sensitiv- 
ity of the relative depth of the blends at 6400A & 6440A with 
effective temperature is very apparent. 



of the mass donor (Fig. g), which clearly shows the rotation- 
ally broadened absorption lines from the donor, especially in 
the wavelength region 6350A- 6540A. Next, this spectrum 
was cross-correlated once again with the velocity standards 
and the relative shifts computed and applied. Then a con- 
stant times each rotationally broadened standard was sub- 
tracted from the average EM Cyg spectrum. The fraction 
chosen was that such that the \ 2 °f the residual spectrum 
was minimized. Plotting the values of \ 2 against the spec- 
tral type of the standard star revealed a minimum at spectral 
types around K3. The circular orbit fit for the K3 standard 
was then shifted out of the contaminant-free EM Cyg spec- 
tra, and the optimal subtraction procedure then repeated to 
check for systematic errors which could appear due to spec- 
tral type mismatch. Initially, the spectral type of the mass- 
donor was estimated by inspe ctio n, in the same manner as 
for the third star (see section 3.1). The absorption line rel- 
ative strengths were noted between wavelengths 6350A and 
6540A as compared to rotationally broadened versions of 
the standard stars. Fig. |^ shows the apparent best position 
to place EM Cyg on the plot, in amongst the standard stars 
whose spectral types range from G8V to M1.5V. Most of the 
M dwarfs for which spectra were obtained have been omit- 
ted from this diagram simply because it became obvious, 
both by inspection and from the results of the initial opti- 
mal subtraction, that they weren't a suitable spectral-type 
match to the mass donor in EM Cyg. The spectral features 
around 6400A and 6440A in particular, were noticeably dif- 
ferent, as was the general shape of the spectrum blueward 
of 6400A. 

After this iterative process had been carried out using 
the standard with the best fit spectral type (broadened ac- 
cording to the measurement of the rotational velocity, which 
is detailed in section B.3), a revised value for the radial ve- 



locity semi-amplitude of the mass donor star was obtained. 
Fig. ^ shows the circular orbit fit to all the data points. The 
parameters of this fit are given in Table [j]. The rms error of 
the fit is lSkms^ 1 . 

Fig. ^| shows the fractional contribution from the mass 
donor to the total light as a function of orbital phase. There 
is no obvious variation in the fraction contributed by the 
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Figure 6. The average Doppler-corrected spectrum of EM Cyg 
after removal of the third star scaled to show the absorption lines 
in the wavelength region of 6390 - 6510A , placed in an appar- 
ently suitable position between spectral types K3V and K5V. The 
sensitivity of the relative depth of the blends at 6400A & 6440A 
with effective temperature is very apparent. 




Figure 7. The 102 measured mass-donor velocities are indicated 
by the data points. The solid line is the fit to these points, details 
of which are given in Table |lj 



mass donor over the orbital period. Also, there appears to 
be no visible decrease around phase 0.5, which indicates 
that irradiation of the donor star by the white dwarf is 
not important. This is confirmed by the lack of any ob- 
vious distortion of the radial velocity curve (see Fig. ^), 
especially around phases 0.25 - 0.75, compared with, for 
example, the distortion seen in the radial velocity curve 
of HS1804+6753 (Billington, Marsh & Dhillon, 1996). A 
weighted least-squares fit to the data points in Fig. ^ gives 
the mean fraction contributed by the mass donor over the 
orbital period to be 0.231±0.005. 



Table 1. Orbital elements and derived parameters 



Element 


Mass DonoiQ 


White Dwar| 




K (kms- 1 ) 


202±3 


170±10 


7 ( kms -1 ) 


-25±2 


-57±7 


T (JDQ) 


2450621.482±0.001 


2443697. 128±0.003 


Msin 3 i (M Q ) 


0.77±0.08 


0.48±0.06 



a Reference: this paper 

6 Reference: Stover et al. (1981) 
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Figure 8. The fractional contribution to the total light by the 
donor star calculated using normalised spectra. For normalised 
spectra, the value obtained from the optimal subtraction process 
for the constant which the standards are multiplied by is the 
fractional contribution of the mass donor star to the total light. 



3.3 Rotational Velocity of the Donor Star 

The rotational broadening of the donor star was estimated 
by shifting the radial velocity fit determined in the previ- 
ous section out of each EM Cyg spectrum, and then co- 
adding the results to produce a spectrum of the star in the 
rest frame of the mass donor (see Fig. ^). Then this aver- 
age spectrum was cross-correlated with the standard stars, 
whose spectra had been artificially broadened by a range 
of velocities (10 - 200 kms -1 in steps of 10 kms -1 ). Both 
the object and standard star spectra were normalised in the 
continuum. The radial-velocity standard spectra were then 
shifted by an amount determined from the correlation. A 
constant times each broadened standard spectrum was then 
subtracted from the object spectrum, to produce a residual 
spectrum. Then a smoothing coefficient was applied to the 
difference spectrum in order to remove any large scale fea- 
tures. The x 2 °f the result was computed in the region con- 
taining the absorption line features (6390 - 6510A ). The re- 
sulting curves of x 2 versus the value of the rotational broad- 
ening used are shown in Fig. for the spectral types from 
Kl to K7. 

From Fig. ^ it is possible to deduce a value for the 
rotational broadening and an estimate of the spectral type 
of the mass donor. Fig. |^ also shows that the systematic 
error due to using different spectral types to obtain values 
for v sini is less than 10 km s - 1 . The minimum x 2 is obtained 
for vsini = 140 ± 3 kms -1 and spectral type K3V. 

Another systematic effect to be taken account when 
calculating the rotational broadening of the mass donor 
is its non-spherical shape. The size of this effect was es- 
timated by computing model absorption profiles including 
the effects of Roche geometry and system parameters ap- 
propriate for EM Cyg. These profiles were then used to 
broaden a K3 standard star. The vsini's of the resulting 
simulated spectra were then estimated using the same pro- 
cedure applied previously to the data. The parameters of 
Table ^| were used to define the Roche lobe (i.e. q = 0.88, 
K\ + K2 = 372 kms -1 ), and the method outlined by Marsh, 
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Table 2. Rotational broadening from model line profiles 



Inclination 


Limb Darkening^ 


v sini 






( km s -1 ) 


69° 





147±4 


69° 


0.5 


140±3 


69° 


1 


130±3 



a Gravity darkening taken to be 0.08 
Mass ratio 0.88 assumed 
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Figure 9. A plot of the x 2 curves obtained using different spec- 
tral type standards, which show that the minima occur at a v sini 
value around 140 kms -1 . 

Robinson & Wood (1994) was used. Table ^ shows the re- 
sults of this process and demonstrates that the distortion 
of Roche lobe geometry has a negligible effect since we re- 
cover vsini — 140 kms -1 for our assumed limb darkening 
of 0.5. However, one can see that increasing the value of the 
limb darkening parameter used to calculate the rotational 
broadening of the lines from the mass donor makes the star 
appear smaller, and introduces a systematic error of the or- 
der 10 kms -1 . In all subsequent analysis the value used for 
the linear limb-darkening parameter was 0.5. 

In addition, the systematic error introduced by alter- 
ing the linear limb-darkening parameter whilst artificially 
broadening the standard stars was also accounted for. Using 
a value of for the linear limb-darkening parameter with 
the best fit standard discovered from Fig. |9[ the measured 
vsini dropped to 135 ± 6 kms -1 . With the limb darkening 
set to 0.5, the value increased to 140 ± 3 kms -1 , an error of 
magnitude 5 kms -1 . 

Accounting for these possible sources of error, we adopt 
a value for the rotational broadening of the mass donor star, 
vsini, of 140 ± 6 kms -1 . From the vsini curves of Fig. bl 
and the relative depths of the absorption lines, the spectral 
type of this star is found to be K3V, and it contributes 23 
per cent of the total light in the Ha wavelength region. 

3.4 Doppler Tomography 

Doppler tomography is an imaging technique which allows 
the distribution of line emission to be mapped in velocity 
space (Marsh & Home, 1988). A thorough analysis of several 
applications of this method to real CV data are given in 
Marsh & Home (1990) and Marsh et al. (1990). 

Fig. [l^ shows the reconstructed Doppler tomogram of 
the Ha emission line in EM Cyg. Any spectra affected by 
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Figure 10. Doppler image of the Ha line. Also plotted are the centre of mass of the donor star at ~V y = 202kms — 1 (upper cross), the 
centre of mass of the system (middle cross), and the centre of mass of the white dwarf (lower cross). The Roche lobe of the mass donor 
star is plotted as well. The dashed line shown is the velocity path of the gas stream. The dash-dotted circle centred on (0,-Ki) is the 
Kcplcrian velocity (430 km s —1 ) at the outer edge of a disc of outer radius O.SRi^ . 



eclipse were not included in the fit. The Doppler image shows 
a diffuse ring-like emission structure, representing the accre- 
tion disc, which if tracking the motion of the white dwarf 
should be centred at position (0, -Ki) in velocity space. 
The inner edge of this ring represents the velocity at the 
outer edge of the disc, which appears to be at a velocity of 
~200 km s" 1 , which is low compared with values of the outer 
disc velocities of other dwarf novae (e.g. Marsh & Home, 
1990); the dashed circle drawn at 430 km s" 1 , is the Keple- 
rian vel ocit y at the outer edge of a disc of radius 0.8R.L! (see 
section 4.2). The apparent 'hole' near the centre-of-mass ve- 



locity of the system is due to the third star. This was not 
removed from the original spectra in order to show its rela- 
tive magnitude and to demonstrate that we can pick up its 
Ha flux. 

The trailed spectra in the left-hand panel of Fig. [Tl| 
shows the set of spectra used to compute the Doppler im- 
age. Visible in Fig. [H] at phase (with the fiducial phase 
being defined as the inferior spectroscopic conjunction of 
the mass donor) is the phenomenon of blue-shifted emis- 
sion being eclipsed before the red-shifted emission (Fig. |l2| ), 
known to be caused by a rotating disc; we use this in the 
next section to constrain the orbital inclination. Although 
there is no obvious S-wave from the mass donor visible in the 
trail, there is emission on the irradiated face of the donor 
visible on the Doppler map (Fig. |l^) . The central panel 
of Fig. [n] shows the data computed from the Ha Doppler 
map. The large region of emission around phase 0.7 appar- 
ent in the actual data, is not completely reconstructed in 



the computed data. Finally, the right-hand panel shows the 
residual formed when the computed data is subtracted from 
the observed data. 



4 DISCUSSION 

4.1 New System Parameters 

The main result of this study has been to discover that the 
spectrum of EM Cyg is contaminated by light from a star 
with a very similar spectral type to its own mass donor. Re- 
moval of the contamination by this star from the spectral 
lines in the region blueward of Ha increases the measured 
semi-amplitude of the radial velocity curve of the mass- 
losing star by over 50kms _1 , effectively removing the prob- 
lem previously associated with EM Cyg: that the system 
should be in a state of dynamically unstable mass transfer. 
Table ^ shows the re-measured orbital elements, and those 
parameters which can be derived from them. 

The range for the inclination i, shown in Table ^| was 
deduced as follows. A simulated Ha trailed spectrum was 
created using the system parameters determined previously 
(namely q and K\ + Ki as 0.88 and 372 km s -1 respectively) 
over the phase range covered by the actual data, i.e. phases 
—0.165 to 0.825. The power law used for the intensity over 
the disc as a function of radius was -1.5, and the velocity 
limits of the plots were ±1500kms _1 . The outer disc radius 
used was 0.8Rl!- The eclipse of the disc in the simulated 
Ha trail becomes visible at inclination angle, i — 61°. The 
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Figure 11. Trailed spectra of the Ha emission line. Emission sites are black on this plot. The left-hand panel shows the observed data 
used to calculate the Doppler tomogram. The middle panel shows the trailed spectrum reconstructed from the tomogram. Finally, the 
right-hand panel displays the residual image formed when the computed data is subtracted from the observed data. 
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Figure 12. Variation of flux from the Ha emission line with orbital phase; 1) The upper plot shows the flux coming from the region 
between 150 and 600 kms -1 2) The lower plot shows the flux coming from the velocity region -150 — 600 kms -1 . The arrows indicate 
the eclipses of the approaching (lower plot) and receding (upper plot) limbs of a rotating disc. 



inclination angle at which the eclipse matches that seen in 
the actual data is i = 67°. At i = 70°, the white dwarf is 
eclipsed, but it is apparent from our data that we are not 
observing an eclipse of the white dwarf (Robinson, 1974). 
This places an upper limit to the inclination of 70°. No part 
of the primary Roche lobe is occulted for inclinations less 
than 61°, thus placing a firm lower limit on i, higher than 
that reached by Robinson (1974). Henceforth, we adopt a 
value of i = 67° ±2° to cover the range where the uncertainty 
should be interpreted as a la estimate. 

The radial velocity semi-amplitude of the white dwarf 
star was checked using the double Gaussian technique de- 
veloped by Schneider & Young (1980). The line wings are 
used to determine the radial velocities because they form 



in an area which is close to the white dwarf and so disrup- 
tion of symmetry by the gas stream and mass donor star are 
minimised. Twin Gaussians (200 kms - FWHM) with sep- 
arations varying from 1200 - 1800 kms -1 were used. This 
procedure gave consistent values of the white dwarf radial 
velocity semi-amplitude K\ = 171 ± 4 kms - , in agreement 
with that obtained by Stover et al. (1981). The systemic 
velocity measured from the radial velocity curve fits was 
7 = —26 ± 5 kms -1 , consistent with that obtained from the 
mass donor analysis. The emission line radial- velocity curves 
were measured to be delayed by 10° ±3° from being exactly 
in anti-phase with with the donor star. 
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Table 3. New orbital and derived system parameters for EM Cygni 



Orbital Element! 



Name 


Value 


Reference 


Period (days) 


0.290909±0.0000042 


Robinson (1974) 


K 2 


202±3kms" 1 


this paper 


Ki 


170±10 kms" 1 


Stover et al. (1981) 


T (JD )(245O621+) 


0.483±0.001 


this paper 


7 ( kms -1 ) 


-25±2kms" 1 


this paper 


a s'mi (cm) 


(1.52 ±0.04) X 10 11 


this paper 


Derived System Parameters 


v sini ( km s ) 


140±6 


this paper 


q (=M 2 /Mi) 


0.88±0.05 


this paper 


M 2 sin 3 i (M ) 


0.77±0.08 


this paper 


Mi sin 3 i (M ) 


0.88±0.05 


this paper 


inclination, i (degrees) 


61° < i < 69° 


this paper 


M 2 (Mq) 


0.99±0.12 


this paper 


Ml (M Q ) 


1.13±0.08 


this paper 


R 2 (R ) 


0.87±0.07 


this paper 



a e = (assumed) 



4.2 Masses 

The measurements of usini and K 2 can be used to deduce 
the mass ratio q = M 2 /Mi. Assuming that the donor star 
rotates synchronously then 



v sini -R 2 (l + i 
K 2 = ~ 



(1) 



The relative size of the donor star is constrained by 
Roche geometry to be 



R.2 



0A9q 



2/3 



0.6<? 2 / 3 + ln(l + q 1 / 3 



for < q < 00, (2) 



Eggleton (1983). Applying these equations we find q = 
0.88 ± 0.05. This is consistent with that found from Ki: 
q = K1/K2 = 0.84 ± 0.06, but we prefer the value deter- 
mined from v sini as Ki measurements often suffer distor- 
tions (Stover, 1981). 

The values obtained here for K 2 , q (derived from the 
rotational broadening measurement) and i can be combined 
with the orbital period obtained by Robinson (1974) to cal- 
culate the masses of each component using: 



Mi 



M 2 = 



PorbK'j(l + q) 2 

27tGsin 3 i 
PorbKl(l + q) 2 q 



(3) 



(4) 



27rGsin 3 i 

The values obtained are shown in Table ^. Assum- 
ing the inclination value concluded in the previous sec- 
tion, i — 67° ±2°, we calculate component masses of 
Mi = 1.12±O.O8M and M 2 = O.99±O.12M . The value for 
the white dwarf mass is consistent with those given for other 
Z Cam systems (Ritter & Kolb, 1998). However, the mass 
for the donor star is close to that of the Sun, which is a G2V 
star, compared with our observed K3 (Allen 1973). Beuer- 
mann et al. (1998) discuss whether the mass donor stars in 
CVs are main-sequence stars, using relationships between 
spectral type and orbital period in a sample of CVs. They 
show that the donor stars in CVs with shorter orbital peri- 



ods fit the 'unevolved' evolutionary tracks well, whilst those 
at longer periods tend to fit 'evolved' models better. EM Cyg 
with its orbital period of 6.98 hours, and mass donor spectral 
type of K3, fits on the plot nearer to the 'moderately evolved 
main-sequence' track, indicating that EM Cyg may have an 
evolved mass donor, because it appears to be oversized and 
over-massive for its spectral type. The radius of the mass 
donor is R2 = 0.87Rq ± 0.07 which compares with O.74R 
for a K5V star, and O.85R for a K0V star (Allen, 1973), 
indicating that it is oversized compared to a main-sequence 
star of the same spectral type. However, the mass deduced 
for the donor star (0.99 ± O.12M ) is not consistent with its 
spectral type according to the current theoretical models of 
Kolb & Baraffe (1999). Their stellar models only become as 
cool as K3 for a donor star mass of £5 O.85M , and this is 
only if the mass loss rate has been rather high (and the mass 
donor is unevolved). In addition, they show that if the mass 
donor is 'evolved' it would have a slightly earlier spectral 
type than a ZAMS star with a mass of 1 M . In summary, 
the new system parameters appear to be inconsistent with 
the current CV mass-donor models. 

We have identified a third star contaminating the spec- 
tral region around Hcv wavelengths in EM Cyg, causing the 
absorption lines in this area to appear narrow, with none of 
the broadening present which is associated with the actual 
mass donor. The spectral type of this third star has been 
identified as K2V - K5V, which is similar to that of the 
mass donor. It also contributes a similar percentage of light 
to the spectrum of EM Cyg as the actual mass donor. The 
third star contributes a fraction 0.160±0.002 to the light of 
EM Cyg at Hcv wavelengths, whilst that contributed by the 
mass donor is 0.231±0.005. Once this third star has been ac- 
counted for, the measurements for the radial velocity semi- 
amplitude increase (from 135±3kms -1 to 202 ±3kms -1 ) 
and the mass ratio of EM Cyg falls below one, removing the 
conflict with the theory. The semi-amplitude of the radial 
velocity curve for the white dwarf is 171 ± 4 kms -1 , which 
agrees with the value 170 ± 10 kms - given by Stover et al. 
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(1981), and which gives a mass ratio of q=0.88±0.05 kms -1 . 
The radial velocity of the centre-of-mass of the binary from 
the absorption line measurement is — 25 ± 2 kms" 1 . The ra- 
dial velocity of the contaminating spectrum is also consistent 
with this value ( — 20±3kms -1 ) which implies that the third 
star and EM Cyg may be physically associated. In addition, 
the similarity of spectral type and fractional contribution to 
the spectrum of EM Cyg indicates that the third star and 
the CV may be at similar distances. However, the field of 
EM Cyg is quite a crowded one, and so a chance superpo- 
sition is also plausible. To address this we acquired some 
direct images on the William Herschel Telescope which have 
a FWHM of 0.7 arcsec, but we still could not resolve the 
third star. The value of asini obtained with the new semi- 
amplitude and mass-ratio values is 1.52 x 10 11 cm. Assuming 
Keplerian motion in the disc, then 

r d _ a GMisin 3 i 
Rlx asini (wdisc sini) 2 ' 

where Kl ± is the distance from the centre of mass of the 
white dwarf to the inner Lagrangian point and is equal 
to 0.51a for q = 0.88. Using the determined values of 
Misin 3 i=O.88±O.O5M , asini = (1.52 ± 0.04) x 10 11 cm, 
and the outer disc velocity, DdiscSini as 305±30kms -1 we 
calculate a value of x^/Rlx = 1.61. This is physically im- 
possible as it implies that the disc must occupy space out- 
side the Roche lobe of the white dwarf. For a more believ- 
able disc radius of 0.8RLJ , the predicted outer disc velocity, 
using equation | is 430 kms -1 . Stover et al. (1981) mea- 
sured the half-separation of the peaks of the H/3 line to be 
390 ± 15 km s - , giving an outer disc radius of 0.98Rli , still 
too large, but not as discrepant as Ha, for which we have 
no explanation. 



5 CONCLUSIONS 

We have found a significant contribution to the spectrum of 
EM Cyg in the Ha wavelength region from a K-type star 
other than the mass donor. When removed from the spec- 
trum, the absorption line radial velocity amplitude increases 
by over 60 kms -1 to 202±3kms -1 bringing the mass ratio 
calculated from the radial velocity semi-amplitudes down to 
0.88 from the value of 1.26 calculated by Stover et al. (1981). 
This solves the 'mystery' surrounding the supposed dynam- 
ical instability EM Cyg would suffer if the mass of the donor 
star was less than O.8M0 and more massive than the white 
dwarf (as deduced by Stover et al., 1981). The radial velocity 
of this 'contaminating' third star (—20 ± 3 kms -1 ) is simi- 
lar to the measured radial velocity of the centre-of-gravity of 
EM Cyg (— 25±2 kms -1 ), which is consistent with the stars 
being physically associated. In addition, the spectral types 
of the two red stars are similar (K2-K5V). The third star 
contributes 16.0±0.2 per cent of the light in this wavelength 
region, 7 per cent less than that contributed by the donor 
star (23.1±0.5 per cent). The rotational broadening of the 
mass donor in EM Cyg is measured to be 140 ± 6 kms -1 , 
and its mass calculated as 0.99 ± O.12M0. Smith & Dhillon 
(1998) plot the mass (in solar mass units) versus the spec- 
tral type of the mass donor in a sample of CVs. EM Cyg 
is on this plot with a spectral type of K5V, and a mass of 
0.76 ± O.IOMq . The new value of M2 puts the system above 



where it should be were it a main-sequence star, confirming 
that the mass donor in EM Cyg must be evolved. If the third 
star is actually physically associated with EM Cyg, then this 
would make it only the second CV known to be in a triple 
system (Reimers, Griffin & Brown 1988). 
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